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The detection of the theoretically expected dark matter is central to particle physics 
cosmology. Current fashionable supersymmetric models provide a natural dark 
matter candidate which is the lightest supersymmetric particle (LSP). Such mod- 
els conmbined with fairly well understood physics like the quark substructure of 
the nucleon and the nuclear form factor and/or the spin response function of the 
nucleus, permit the evaluation of the event rate for LSP-nucleus elastic scattering. 
The thus obtained event rates are, however, very low or even undetectable. So it is 
imperative to exploit the modulation effect, i.e. the dependence of the event irate 
on the earth's annual motion. In this paper we study such a modulation effect both 
in non directional and directional experiments. We calculate both the differential 
and the total rates using symmetric as well as asymmetric velocity distributions. 
We find that in the symmetric case the modulation amplitude is small, less than 
0.07. There exist, however, regions of the phase space and experimental conditions 
such that the effect can become larger. The inclusion of asymmetry, with a realistic 
enhanced velocity dispersion in the galactocentric direction, yields the bonus of an 
enhanced modulation effect, with an amplitude which for certain parameters can 
become as large as 0.46 



1 Introduction 

It is known that that dark matter is needed to close the Universe 0,i. It is also 
known that one needs two kinds of dark matter. One composed of particles 
which were relativistic at the time of structure formation. These constitute 
the hot dark matter component (HDM). The other is made up of particles 
which were non-relativistic at the time of freeze out. This is the cold dark 
matter component (CDM). The COBE data 1 suggest that CDM is at least 
60% □. Druthe other hand recent data from the Supernova Cosmology Project 
suggest u ,Q that there is no need for HDM and the situation can be adequately 
described by 51 < 1, e.g. ^cdm = 0.3 and JIa — 0.6. In a more recent analysis 
TurnerQ gives — 0.4. 

Since the non exotic component cannot exceed 40% of the CDM there 
is room for the exotic WIMP's (Interacting Massive Particles). Recently the 
DAMA experiment □ has claimed the observation of one signal in direct detec- 
tion of a WIMP, which mifh better statistics has subsequently been interpreted 
as a modulation signal ll^. 
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In the currently favored supersymmetric extensions of the standard model 
the most natural WIMP candidate is the LSP, i.e. the lightest supersymmetric 
particle, whose nature can be described in most supersymmetric (SUSY) mod- 
els to be a Majorana fermion, aJinmr combination of the neutral components 
of the gauginos and Higgsinos 

Since this particle is expected to be very massive, > SOGeV, and 
extremely non relaliviatic with average kinetic energy T < lOOKeV, it can be 
directly detected only via the recoiling of a nucleus (A,Z) in the elastic 

scattering process: 

X + {A,Z) ^ X + {A,zr (1) 

(x denotes the LSP). In order to compute the event rate one proceeds with 
the following steps: 

1) Write down the effective Lagrangian at the elementary particle (quar 
level obtained-in the_itamework of supersymmetry as described in Refs 
Bottino et alM andN. 

2) Go from the quark to the nucleon level using an appropriate quark model 
for the nucleon. Special attention in this step is paid to the scalar couplings, 
which dominate the coherent part of the cross section and the isoscalar axial 
jTii.rr/jnt, which, as we will see, strongly depend on the assumed quark model 



3) Compute the relevant nuclear matrix elementsEjTJ using as reliable as 
possible many body nuclear wave functions hoping that, by putting as accurate 
nuclear physics input as possible, one will be able to constrain the SUSY 
parameters as much as possible. 

4) Calculate the modulation of the cross sections due to the earth's re¥Aj 
lution around the sun by a folding procedure assuming some distribution BEa 
of velocities for LSP. 

The purpose of our present review is to focus on the_last point of our 
above list along the lines suggested by our recent letter E3, expanding our 
previous results and giving some of the missing calculational details. For the 
reader's convenience, however, we will give a brief description on the basic 
ingredients on how to calculate LSP-nucleus scattering cross section, without 
elaborating on how one gets the needed parameters from supersymmetry. For 
the calculation of these parameters from representative input in the restricted 
SUSY parameter space, we refer the reader to the literature, e.g. Bottino et al. 
l3, Kane et al. , Castano et al. and Arnowitt et al. Ej. Then we will specialize 
our study Jn the case of the nucleus ^^^/ which is one of the most popular 
targets E£l LI. To this end we will include the effect of the nuclear fomi factors. 
We will consider both a symmetric Maxwell-Boltzmann distcjbutionEI as well as 
asymmetric distributions like the one suggested by DrukierllS. We will examine 
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the effect modulation in the directional as well as the non directional detection, 
both in the differential as well as the total event rates. We will present our 
results a function of the LSP mass, m^, for various detector energy thresholds, 
in a way which can be easily understood by the experimentalists. 



2 The Basic Ingredients for LSP Nucleus Scattering 

Because of lack of space we are not going to elaborate here further on the 
construction of the effective LaKnangian .derived from supersymmetry, but re- 
fer the reader to the literature til't30l23S The effective Lagrangian can be 
obtained in first order via Higgs exchange, s-quark exchange and Z-exchange. 
We will use a formalism which is familiar from the theory of weak interactions, 
i.e. 

Leff = -^{ixil^l5Xi)Jx + iXiXi)J} (2) 

where 

Jx = Njx{f^ + fl^r3 + fh5 + fh5T3)N (3) 

and 

J = N{f° + flT,)N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, Xil^Xi = (identically). 
The parameters fy, fy, f\, fg: fs depend on the SUSY model employed. 
In SUSY models derived from minimal SUGRA the allowed parameter space 
is characterized at the GUT scale by five parameters, two universal mass pa- 
rameters, one for the scalars, ttiq, and one for the fermions, mi p, as well as 
the parameters tanf3, one of Aq and "''^ and the sign of /i o. Deviations 
£cpm universality at the GUT scale have also been considered and found useful 
Ej. We will not elaborate further on this point since the above pararnetj?r.s in- 
volving universal masses have already been computed in some modelsFHri and 
jeffects resulting from deviations frnm universality will be published elsewhere 
Eil (see also Arnowitt et al in Ref. ES and Bottino et al in Ref. E2l). For some 
choices in the alloss/ed parameter space the obtained couplings can be found in 
a previous paperEj. i_. 

The invariant amplitude in the case of non-relativistic LSP can be cast llil 
in the form 



, ^ EfE, - ml + p. . p, ^ ^ 
mi 

c^/32|Jo|2 + |J|2 + |j|2 (5) 



3 



where rrix is the LSP mass, | Jo| and |J| indicate the matrix elements of the 
time and space components of the current Jx of Eq. (H), respectively, and J 
represents the matrix element of the scalar current J of Eq. (||). Notice that 
I Jop is multiplied by 0^ (the suppression due to the Majorana nature of LSP 
mentioned above). It is straightforward to show that 



J^^A^\F{ce)? (7) 

with F(q^) the nuclear form factor and 

A A 
no(q) - ^a(J)e-^'^-^ ^M) = E ^(j>3(j>-*'' ''^ (9) 

where crO), T3(j), x^- are the spin, third component of isospin (rajp) = \p)) 
and coordinate of the j-th nucleon and q is the momentum transferred to the 
nucleus. i_. 
The differential cross section in the laboratory frame takes the formEJ 



% = ^(^) wi JoHi - ^^^e] + |jp + (10) 

a\l IT niN (1 + ^7) 



where m^r is the proton mass, rj = mx/m^A, ^ = Pi • q > (forward 
scattering) and 

<Jo = i^iGprnNf ^0.77 X IQ-^^ cm^ (11) 

ZTT 

The reduced mass /i^ is given by 

Mr = — ^ 12 
1+7] 

For the evaluation of the differential rate, which is the main subject of the 
present work, it will be more convenient to use the variables {v,u) instead of 
the variables (w,0- Thus integrating the differential cross section, Eq. (p^, 
with respect to the azimuthal angle we obtain 

'^^("' ^) = 7:7^T-T2 [(^S +^V^ F^{u)) + tsp^nFll\ (13) 
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with 

ts = a,{J^r{A'[{fs~fs^^f] (14) 
E,,,„ = aoi^r [flMO))'^, + 2/l/Af^o(0)17i(0)|4l + if\^m)' ] 

niM J:'ll(u) ril{Uj 

(15) 

We should remark that even though the quantity Espm can be a function of u, 
in actual practice it is indepenent of u. The same is true of te less important 
term Ey In the above expressions F{u) is the nuclear form factor and 

are the spin form factors with 

u = q%^/2 (18) 

b being the harmonic oscillator size parameter and q the momentum trasfer to 
the nucleus. The quantity u is also related to the experimentally measurable 
energy transfer Q via the relations 

The detection rate for a particle with velocity v and a target with mass m 
detecting in the direction e will be denoted by R{—^ e). Then one defines the 
undirectiona rate Rundir via the equations via the equations 

Rundir = -JT ^ A ^-^x I + I '"■^y I + I "^-e^ IJ (20) 

at ArriN 

p(0) = 0.3GeV/cm^ is the LSP density in our vicinity. This density has to be 
consistent with the LSP velocity distribution (see next section). 
The differential undirectional rate can be written as 

dRundir = da{u,v)[\ V.Cx I + I v.Sy I + I V.Ez \] (21) 

where da{u, v) is given by Eq. ( |l^) 
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The directional rate in the direction e takes the form: 

Ra,r - e) - - e) = ^ -^„.e a{u, v) (22) 

ArriN 

and the corresponding differential rate is given by 

dRdir = v.e da{u,v) (23) 

niy, Attin 

3 Convolution of the Event Rate 

We have seen that the event rate for LSP-nucleus scattering depends on the 
relative LSP-target velocity. In this section we will examine the consequences of 
the earth's revolution around the sun (the effect of its rotation around its axis is 
expected to be negligible) i.e. the modulation effect. This can be accomplished 
by convoluting the rate with the LSP velocity jxlistribution. Hitherto such a 
consistent choice can be a Maxwell distribution □ 

/(«') - iV^vor^e-^^'/^o')" (24) 



Wo = V{2/3)y) = 220Km/s (25) 

i.e. vo is the velocity of the sun around the center jof the galaxy. In the present 
paper following the work of Drukier, see Ref. ES, we will assume that the 
velocity distribution is only axially symmetric, i.e. of the form 



with 



fi{v , A) = exp[{ ^ ] (27) 

/2 w , Vesc, A) = exp[ "-^ (28 

where fesqjf the escape velocity in the gravitational field of the galaxy, Vesc ~ 
625Km/ si^. In the above expressions A is a parameter, which describes the 
asymmetry and takes values between and 1 and N is a proper normalization 
constant given by 

1 - 1 [e./(y...)-e-(^+^)e^^^^(^^^-=) 



N{X,yesc) X + V i^A 

Q-y^sa 2 ), „2 erf {i \/X Vesc) ^ 

A y/TT i VA 



with Uesc — and erf(x) the error function given by 

er/(x) = ^ / die"*' (30) 
V'^ Jq 

For j/esc ^ 00 we get the simple expression N^^ A + 1 

The z-axis is chosen in the direction of the disc's rotation, i.e. in the 
direction of the motion of the the sun, the y-axis is perpendicular to the plane 
of the galaxy and the x-axis is in the radial direction. In the above frame we 
find that the position of the axis of the ecliptic is determined by the angle 
7 « 29.80 (galactic latitude) and th^^azimuthal angle lu = 186.3" measured on 
the galactic plane from the z axis 113. Thus, the axis of the ecliptic lies very 
close to the y, z plane and the velocity of the earth around the sun is 

'Ve = + v\ = vq + vi{ sinaSi. ~ cosacosjy + cosa sin"f z) (31) 

where a is the phase of the earth's orbital motion, a = 2Tr(t — ti)/TE, where 
ti is around second of June and Te — lyear. 

We are now in a position to express the above distribution in the laboratory 
frame, i.e. 

f{v, X, Ve) = h{v, VE, A) - f4,{v, Vesc, A) (32) 

with 



(1 + X){{vy + vicosjsina)^ + [v^ + wq + visin^cosaf') 



h{v,VE,X) = exp[- 
X ea;p[— 



. , , , r + X((vy + vicosjsina)'^ + (v^ + ^-^o + visin^fcosa)'^) , 
/4(w, Vesc, A = exp[ — ■ '-^ —] 

(34) 



4 Expressions for the Diflferential Event Rate in the Presence of 
Velocity Dispersion 



We will begin with the undirectional rate. 
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Expressions for the Undirectional Differential Event Rate 
The mean value of the undirectional event rate of Eq. ( p^ ) , is given by 

/dRundir\ p{0) m f -1,1 -1,1 - i-idcr{u,v) . 
{ — ^ ) = 1 / f(v,VE)[\ v.e^ I + I v.ey I + I v.e^ \] — d'v 

(35) 

From now on we will omit the subscript undir in the case of the undirectional 
rate. The above expression can be more conveniently written as 

^dR.piOln^^^d^ (36) 
\ du / Am n du 

where 

.dE f [\\v.e4 + ||i;.e^|| + \\v.e,\\] ^^^^ 3 

{-!-/ = / Tf^f J(v,ve) — a V (67) 

du J Viv ) du 

It is convenient to work in spherical coordinates. But even then the angular 
sun is small. Thus introducing the parameter 

(5 = — = 0.27, (38) 
Wo 

expanding in powers of S and keeping terms up to linear in it we can manage 
to perform the (f> integration using standard contour integral techniques and 
express the result in terms of the two modified Bessel functions Imi^^i^^f^)) 
with t = cos9 and m=0,l- Thus the angular integration of Eq. ^ yields 

X exp[-{y^ + 1)(1 + A)]A,(A, y) - exp[-{yl, + Xy^)]A-{X, y), i = 1, 2 

(39) 

where A^, Aj come from /a, fi respectively. We find 

Ai(A,y) = Ai,i(A,y)+Ai,2(A,y)+Ai,3(A,y) (40) 

with ^ 

Ai,3(A, y)= J ^ dtexp[iCV2 - 2(A + l)yt)] \ t \ /o(CV2) (41) 

Ai,2(A,y) = ^ f\texp[{^2{X + l)yt)]{l-tY/'^^:I^ (42) 

vi" J-i « C 



Ai,i(A,y) = ^ dtexp[{^2{\ + l)yt)]il - ^2)i/2 e^P(C^) er/(C) ^^g^ 

where C = y (A (1 — t^))^/^ and the function erf[iC,) given by Eq. 
Furthermore 

A2(A,y) = -(A + l)A;'(A,y) (44) 

where Aj^ (A, y) is obtained from Ai(A, y) by adding in the integrands the extra 
factor ty+1. We should mention that in the last integral we have ommitted the 
numerical factor Scosasin'^ . Note that in the case A2(A, y), which is asociated 
with the modulation amplitude, only the z-component of the velocity in the 
exponential contributes. Hence the dependence on the earth's phase is cosa. 
The formulas for the second term in Eq. ^ for Aj(A,y) are obtained by a 
mere replacement of the expression A + 1 by A. In all the above expressions 
y — {v/vq) (not to be confused with the y-coordinate) . 

It is convenient to separate out the asymmetric contribution from the usual 
one by writing 

2y2Ai(A, y) = #o(A, (A + l)2y) + Go(A, y) (45) 

2y2A2(A,2/) = A(A,(A + l)2j/) + Gi(A,2/) (46) 

The functions Fi have been obtained by considering the leading non vanishing 
term in the zeta expansion of the integrands of the expressions (^l|)-(^3|). Thus 

Go(0,y) = , Gi(0,2/) = (47) 

Fq{\x) = (A + iy^[x sinh{x) - cos{x) + l + x Ii{x)] (48) 



Fi{\,x) = (1 + A)"2 [(2 + A)((xV(2(2 + A)) + 1) cosh{x) - x sinh{x) - 1) 
+ x'' l2[x\-{\ + l) X h{x)] (49) 

note that here a; = (A + 1)2?/. Im{x) is the modified bessel function of order 
m. The funcions G cannot bo obtained analytically, but they can easily be 
expressed as a rapidly convergent series in y = which will not be given 
here. 

Similarly 

G;(A,y) =2y2A;(A,y) , i = l,2 (50) 
Thus the folded non-directional event rate takes the form 

= tsFo (u) + ^tvFi (u) + tsp^nFsp^n {u) (51) 
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where the — S,V, spin are given by Eqs. (|lj)- (p^). 

The quantities Fq, Fi, Fspin are obtained from the corresponding form fac- 
tors via the equations 

Fk{u) ^ F\u)^k{u) ^\l^^f , fc-0,1 (52) 
2fc + 1 

FspM^Fix{u)-^o{u)a^ (53) 

(54) 

§fc(w) = [-0(0), fc («%/") + 0.135 cos Q;'(/'(i)^fe(ayu)] (55) 



with 



and 



1 



(56) 



^(^i),k{^) = N{yesc:\)e ^(e - exp[~yU^[i^^ ,^{x)) (57) 
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<i>(,).fe(a;) = -= / dyy^'-'expi-il + \)y')){Fi{\ (A + l)2y) + Gi{\ y))) 

VOTT Jx 

I-;,) ,(x) = -|=^''°'dyy2fc-ie^p(„Ay2))G;(A,2/)) (59) 
The undirectional differential rate takes the form 

(^) = RtT{um + cosai/(M))] (60) 

In the above expressions R is the rate obtained in the conventional approach 
IliJ by neglecting the folding with the LSP velocity and the momentum transfer 
dependence of the differential cross section, i.e. by 

p(0) m 



+ (61) 

where = S,V, spin have been defined above, see Eqs (|lj) - ([T6|). 

The factor T{u) takes care of the u-dependence of the unmodulated dif- 
ferential rate. It is defined so that 

duT{u) = 1. (62) 
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i.e. it is the relative difFrcntial rate. Umin is determined by the energy cutoff 
due to the performance of the detector. Umax is determined by the escape 
velocity Vesc via the relations: 

Vesc ff\'i\ 

Umax — ^ V^*^/ 

On the other hand H[u) gives the energy tranfer dependent modulation ampli- 
tude. The quantity t takes care of the modification of the total rate due to the 
nuclear form factor and the folding with the LSP velocity distribution. Since 
the functions Fq{u),Fi and Fspm in principle have a different dependence on 
u, the funcions T(u), H[u) and t in priciple depend on the SUSY parameters. 
If, however, we ignore the small vector contribution and assume (i) the scalar 
and axial (spin) dependence on u is the same or (ii) only one mechanism (S, 
V, spin) dominates the, parameter R contains the dependence on all SUSY 
parameters. The other factors depend only on the LSP mass and the nuclear 
parameters. More specifically considering only the scalar interaction we get 
R Rg and 

tT{u)^a^F\u)i^^o)AaV^) (64) 
For the spin interaction we get a similar expression except that R — s- Rspin 
and F^ — s- -Fp.p' . Finally for completeness we will consider the less important 
vector contribution. We get R Ry and 

1 2?/ + 1 ~ ^ 2a2 

The quantity T{u) depends on nuclear physics through the form factors or the 
spin response functions and the parameter a. The modulation amplitude takes 
the form 

,/(,)^0.135tw^^,/ = l,3 (66) 
V'(o),/c(a\/w) 

Thus in this case thj2. H[u) depends only on a^/u, which coincides with the 
parameter x of Ref. E3, i.e. only on the momentum transfer, the reduced mass 
and the size of the nucleus. 

Returning to the differential rate it is sometimes convenient to use the 
quantity T{u)H{u) rather than iJ, since H(u) may appear artificially increasing 
function of u due to the decrease of T{u) ( in obtaining H(u) we have divided 
T(u)) 

Before concluding this subsection we should mention that the above an- 
gular integrations can also be done even if the velocity distribution is triaxial. 
We will not explore this further since one has too many parameters. 



t Tiu) = F'^{u)[ij(Q)^i{a^/u) - 7^777^ 77:^-;TjM?/'(o),o(aVw)]^ (65) 
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4-2 Expressions for the Directional Differential Event Rate 

The mean value of the directional differential event rate of Eq. (|23|), is defined 
by 

/dR\ p(0) m f ^, . daiu.v) ,n 

{ — ) / f{v,VE)v.e \ ' dPv 67 

\ au / dir Am,]^ J au 

where e is the unit vector in the direction of observation. It can be more 
conveniently expressed as 

(f) -^i^^if) (68) 
\ au / dir AmM au dir 

where 

dE f v.e da{u,v)^ 

{-r-)d^r^ / ^— f{v,VE) — d-'v (69) 

du J y'(v^) du 

Working as in the previous subsection, i.e by expanding in powers of 5 and 
keeping terms up to linear in it we can manage to perform the (p integration 
using standard contour integral techniques and express the result in terms of 
the two modified Bessel functions Im{^^{^ ~ ^^)) with t — cos6 and m=0,l. 
Thus the angular integration of Eq. ^ yields 

M,(A,y) = 27r 

X exp[-{y'' + 1)(1 + A)]A,(A, y) - exp[-{yl, + Xy^)]A-i\, y),i^ 1,4 

(70) 

where A^, A — i come from /a, f^ respectively and are given by 
Ai(A,j/)= dtexp[-iiX/2)y^il-t^) + 2iX+l)yt)]tIo{{X/2)y^il-t^)) (71) 



A2(A,2/)= / dtexp[^{{X/2)y\l-t') + 2iX+l)yt)]t'loiiX/2)y'il^t')) (72) 



A3(A,y) = J dtexp[-iiX/2)y'{l-t') + 2{X+l)yt)]{l-t')Io{{X/2)y'{l^t')) 

(73) 

A,{X,y) = J dtexp[-{{X/2)y^{l-t^) + 2{X+l)yt)]{l^t')hi{X/2)y^l-t^)) 

(74) 

and analogous expressions for A with the mere replacement of the expression 
A + 1 by A. Again in the above expressions y=(u/fo) (not to be confused 
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with the y-coordinatc). The above integrals can be expressed in terms of 
hypergeomctric functions as follows 

, , ^ 2(2(A+ ,1 2fc + 5 , 



fe=0 



UX, y) = iX + l)Xf (2fe + l)(2fc + 3)(2fc + 5)((2fc)!) ^^ (2' ^ 



fc=0 



, ^ 2(2Ay)(2fc+i) 1 2fc + 5 ^ , 

= - (2fc + 3)((2fc + l)!) ^^ ^2 ' -^'^y 



k=0 



(78) 
(79) 



, 2(2Ay)(2fe) .1 2fc + 5 , , 

^) = ^ ^ (2feV3)((2fc)!) ^^ (2' 



fc=0 



) (80) 



,'0,V_xV 2(2A2;)(^^) ,1 2fc + 5 

A3(A,yJ - A 2^ ^2^, ^ ^^^2^ ^ 3)((2/fc)!) 2 4' ' ' 2 



/c=0 



.As') (81) 



'^•(^•'> - '^^)'g p.+ l)(2'';^|^'+5)((2m) ^4'^-^'''' <«^) 

It is more convenient to define the functions Fi ang G^, i = 0,4 as follows 

~2y^Ai{X,y)^Fo{2iX+l)y) + Go{X,y) (83) 

2y2(Ai(A, y) + yA2(A, y)) = Fi(A, 2(A + 1)) + Gi(A, y) (84) 

4y3(A3(A, y) - A4(A, y)) = F2(2(A + 1)) + G2(A, y) (85) 

4y3(A3(A, y) + A4(A, y)) = F3(2(A + 1)) + G3(A, y) (86) 

The functions Fi are obtained by keeping the leading terms in the expansion 
of the confluent hypergeommetric functions of Eqs. (|75|)- (78). Thus we find 

Ftix) = Xcoshx - sinhx , i = 0, 2, 3 (87) 
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Fi(A,x) = 2(1 -A) [{^-^±^ + l).sinhx-Xco.shx 

The purely asymmetric quantities Gi satisfy 

G,(0,y) = 0, (89) 

They are expressed as rapidly convergent series in y, but they are not going to 
be given here. Similarly we define the functions G via the equations 

G'a{X,y) = -2i/A[{X,y) (90) 

G[ (A, y) = 2y\A[{X, y) + y^^{\ y)) (91) 

G2{\,y) = 4y3(A;(A,y) - A'4(A,y)) (92) 

G;(A,y) = (93) 
Thus the folded directional event rate takes the form 

where the , i = y, spin are given by Eqs. (H)- (|l|). The quantities 

Fq, Fi, Fgpin are now obtained from the corresponding form factors via the 
equations 

Fk{u)^F^u)^,{u)^l±^,k = 0,1 (95) 
Zk + 1 

Fsp^n{u) = Fu{u)^oiu)a^ (96) 

= ^[(V'(o).fe(a\/u) + 0.135 cosQ;?/;(i)^fc(aVw))ez.e 

— 0. 117 cos aV'(2), fc (av^)ey. e + 0.135 sin a'(/'(3)_fe(av^)ea;.e] (97) 

with 

i^(i).k{x) = A^(2/esc, A)e-^(e-i$(i),fc(a;) - exp[-y2j*'(0,fc(2^)) (98) 
$(;),,(x) = ^ r^'^ dyy'''-'expi-il + X)y^))iFii2y) + Gii\,y))) (99) 

V DTT 

<i>;;),,(a;) = ^ J'"^ dyy^''-'expi~Xy^))G'iiX,y)) (100) 
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If we consider each mode (scalar, spin vector) separately the directional rate 
takes the form 



dR R 

{-j-)dir ~ "2" R^\{1 + COS aHi{u))ez.e — COS aH2{u)ey.e + sin aH3[u)ex-e\ 

(101) 

In other words the non directional differential modulated amplitude is de- 
scribed in terms of the three parameters, Hi{u), 1=1,2 and 3. The unmodulated 
one is R^{u), which is again normalized to unity. It is the relative differen- 
tial rate, i.e. the differential rate divided by the total rate, in the absence of 



modulation, i.e. The parameter i° entering Eq. (101) takes care of whatever 
modifications are needed due to the convolution of the non modulated total 
rate with the LSP velocity distribution in the presence of the nuclear form 
factors. 



From Eqs. - ( 101 ) we see that if we consider each mode separately 



the differential modulation amplitudes H(l) take the form 

ff;(w) = 0.135— ^y- — — , / = 1,3 ; i72(w) = 0.117-^y- — — (102) 

Thus in this caseJJjie Hi depend only on ay^, which coincides with the pa- 
rameter X of Ref. E3. We note that in the case A = we have H2 = 0.117 and 
H3 = 0.135, so that there remains This means that, if we neglect the coherent 
vector contribution, which, as we have mentioned, is small. Hi essentially de- 
pends only on the momentum transfer, the reduced mass and the size of the 
nucleus. 

Returning to the differential rate it is sometimes convenient to use the 
quantity Ri rather than Hi defined by 

Ri = R" Hi, I = 1,2,3. (103) 

The reason is that Hi, being the ratio of two quantities, may appear superfi- 
cially large due to the denominator becoming small, once again if one mech- 
anism dominates the parameters Rq and Ri are independent of the particular 
SUSY model considered, except the LSP mass. In fact we find for the scalar 
interaction we get R ^ Rs and 

t° R°{u) = a^F\u)^l;l^\ay^) (104) 

For the spin interaction we get a similar expression except that R — > Rspin 
and -F^ Fp^p' . Finally for completeness we will consider the less important 
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vector contribution. We get R Rv and 

R\u) = F\u)[i,f\aV^) - ^^^^^«vr(aV^)]^ (105) 

The quantity i?o depends on nuclear physics through the form factors or the 
spin response f unct ions. 

Equation ( 101 ) deviates from the simple trigonometric expression. a The 
dependence on the phase of the earth is complicated. If we imagine, however, 
that one can sum up all three directonal rates, with the inclusion of H2 and 
H^, the maximum does not occur at a = 0, but at a = an with 

o^H^tan-^ ^'^""l ] (106) 

and the modulation at this value of the phase of the earth takes the value 

H,nax = [{Hi + i?2)' + Hi] i (107) 

There exists one minimum at a = tt, i.e. around Dec. 2 and takes the value 

Hmin — H2 — Hi (108) 

Whenever Hi > H2 there exist two more minima a — 7r/2 and 3(7r/2), equal 
to and two secondary maxima. In all cases considered in this work > 



H2 — Hi so that the interesting quantities are given by Eqs(107 - |10q). In 



any case it is useful to know the difference between the maximum and the 
minimum, which takes the form 

Hra = [(i^i + H2f + Hl]^ - Min{Hi - H2, H3) (109) 
5 The Total Modulated Event Rates 

Once again we will distinguish two possib ilitie s, namely the directional and 
the non directional case. Integrating Eq. (101) we obtain for the total undi- 
rectional rate 

R^ Rt[{l + h{a,Q^„,)cosa)] (110) 

where Qmin is the energy transfer cutoff imposed by the detector. The modu- 
lation of the non-directional total event rate can be described in terms of the 
parameter h. 

The effect of folding with LSP velocity on the total rate is taken into 
account via the quantity t. The SUSY parameters have been absorbed in 
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R. From our discussion in the case of differential rate it is clear that strictly 
speaking the quantities t and h also depend on the SUSY parameters. They 
do not depend on them, however, if one considers the scalar, spin etc. modes 
separately. 

Let us now examine the directional rate. Integrating Eq. ( p5|) we obtain 

Rdir = R (i"/2) [(1 + hi{a, Qmin)cosa)e ^.e 

- /i2(a, Qrnin) cosae y.e + hsia, Qmin) sinae ^-e] (111) 

Furthermore if we somehow manage to measure the directional rate in all 
directions we obtain: 

Rdir,aU = R (^°/2) [l+/ii(a, Qmm) co3a+h2{a, Q,nin) \cosa\+hz[a, QnUn) \sina\] 

(112) 

We see that the modulation of the directional total event rate can be described 
in terms of three parameters hi 1=1,2,3. In the special case of A = we 
essentially have one parameter, namely hi, since then we have /12 = 0.117 and 
hs = 0.135. 

The effect of folding with LSP velocity on the total rate is taken into 
account via the quantity All other SUSY parameters have been absorbed 
in R, under the same asumptions discussed above in the case of undirectional 
rates. 

Given the functions/i; (a, Qmm) one can plot the the expression in Eq. 
(112) as a function of the phase of the earth a. For a gross description one can 
follow the procedure outlined above making the substitution H h. Thus the 
maximum occurs at ±ah with 

The difference between the the maximum and the minimum is now given by 

/i„ = [{hi + h2f + hl]^- Min{hi - h2,h3)\ (114) 
In all cases considered here h^ > \hi — h2\ 



6 Discussion of the Results 

We have calculated the differential as well as the total event rates (directional 
and non directional) for elastic LSP-nucleus scattering for the target 
including realistic form factors. Only the coherent mode due to the scalar 
interaction was considered. The spin contribution will appear elsewhere. Spe- 
cial attention was paid to the modulation effect due to the annual motion of 
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the earth. To this end we included not only the component of the earth's 
velocity in the direction of the sun's motion, as it has been done so far, but 
all of its components. In addition both spherically symmetrica as well only 
axially symmetricEa LSP velocity distributions were examined. Furthermor we 
considered the effects of the detector energy cutoffs, by studying two typical 
cases Qmin = 10 and 20 KeV both on the modulated and the unmodulated 
aplitudes. We focused our attention on those aspects which do not depend on 
the parameters of supersymmetry other then the LSP mass. 

The parameter R, normally calculated in SUSY theories, was nnt consid- 
ered in this work. The interested reader isjjejferred to the literature Ej and, 
in our notation, to our previous workll^l ' tHi E3. 



6.1 The Undirectional Rates 

Let us begin with the total rates, i.e. the quantities t and h. In Table 1 we 
show the dependence of h on the components of the earth's velocity, in the 
symmetric case(A = 0). We see that the modulation amplitude increases for 
about 50 % when all components of the earth's velocity are included. 

In Table 2 we show how the quantities t and h depend on the detector 
energy cutoff and the LSP mass for the symmetric case. In tables 3 and 4 we 
show the same quantities for A — 0.5 and A = 1.0 respectively. From these 
tables we see a dramatic incease of the modulation when the realistic axially 
symmetric velocity distribution is turned on. This means that the modulation 
amplitude can be exploited by the experimentalists. We further notice that the 
modulation amplitude increases somewhat with cutoff enery. This is due to 
the fact that the modulation amplitude decreases less rapidly with the cutoff 
energy Qmin than the unmodulated amplitude. This effect may be of use to 
the experimentalists, even though it occurs at the expense of the total rate. 



Let us now examine the differential rates, which are described by the func- 
tions T{u), H{u) and T{u) x H{u). These are shown for various LSP masses 
and Qmtn in Fig. 1 (A = 0.0), Fig. 2 (A = 0.5) and Fig. 3 (A = 1.0) We 
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Table 1: The dependence of the moduahion amplitude h on the velocity of the earth in the 
symmetric case A = 0. and Qmin = 0. 









LSP 


mass 


in GeV 






Velocity 


10 


30 


50 


80 


100 


125 


250 


z-com 
all 


0.0453 
0.0723 


0.0320 
0.0558 


0.0179 
0.0383 


0.0075 
0.0252 


0.0041 
0.0208 


0.0015 
0.0173 


-0.0033 
0.0112 



Table 2: The quantities t and h for A = in the case of the target 53/^ for various 
LSP masses and Q,nin in KeV (for definitions see text). Only the scalar contribution is 
considered. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t 


0.0 


1.599 


1.134 


0.765 


0.491 


0.399 


0.328 


0.198 


h 


0.0 


0.072 


0.056 


0.038 


0.025 


0.021 


0.017 


0.011 


t 


10. 


0.000 


0.276 


0.307 


0.236 


0.200 


0.170 


0.108 


h 


10. 


0.000 


0.055 


0.028 


0.014 


0.010 


0.007 


0.001 


t 


20. 


0.000 


0.058 


0.117 


0.110 


0.098 


0.086 


0.058 


h 


20. 


0.000 


0.084 


0.044 


0.024 


0.017 


0.013 


0.005 



remind the reader that the dimensionless quantity u is related to the energy 
transfer Q via Eq. ( pO[ ) with Qo = 60KeV for ^^^/. The curves shown corre- 
spond to LSP masses as follows: i) Solid line = 30 GeV. ii) Dotted 
line -<=4' — 50 GeV. iii) Dashed line = 80 GeV. v) Intermediate 
dashed line ^=> = 100 GeV. vi) Fine solid line <^=^ — 125 GeV. vi) 
Long dashed line 4=^ — 250 GeV. If some curves of the above list seem 
to have been omitted, it is understood that they fall on top of vi). Note that, 
due to our normalization of T, the area under the corresponding curve is unity. 
This normalzation was adopted to bring the various graphs on scale since the 
absolute values may change much faster as a function of the LSP mass. 
In order to understand the dependence of the total and differential rates on A, 
we will examine the functions f{y), which are equal to the quantity A^(A, y) 
multiplied by the integrand of Eq. (|5^) . THe latter crucially depends on the 
functions ^^(A, 2(A + l)y) and Gi{X,y) , i — 1,2. These functions f{y) are 
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T( 
1.2 

1 

0.8 

0.6 
. 4 
0.2 



u) 



0.1 0.2 0.3 0.4 



T(u) *H(u) 

0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 



0.1 0.2 0.3 0.4 O.B 



(A = 0.0, Qmin = 00.) 



(A = 0.0, Qmin = 00.) 





0.2 0.3 0.4 



0.5 0.6 



0.2 0.3 0.4 0.5 0.6 



(A = 0.0, Qmin = 10.) 



(A = 0.0, Qmin = 10.) 



T(u) 





0.4 0.5 0.6 0.7 



0.4 0.5 0.6 0.7 



(A = 0.0, Qmin = 20.0) 



(A = 0.0, Qmin = 20.) 




0.1 0.2 0.3 0.4 O.S 



(A = 0.0, independent of Qmin) 



Figure 1: The quantities T, H and TH entering the undirectional differential rate for A = 0.0 
and various values of energy cit off in GeV . For definitions see text. The energy trasfer Q 
is given by Q = uQo, Qo = GOKeV 



T( 
2 

1.5 

1 

0.5 



u) 



0.1 0.2 0.3 0.4 O.S ^ 



T(u) ' 

0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



H(u) 



0.1 0.2 0.3 0.4 O.B 



(A = 0.5, Qmin = 00.) 



(A = 0.5, Qmin = 00.) 




*H(U) 




0.2 0.3 0.4 0.5 0.6 



0.2 0.3 0.4 0.5 0.6 



(A = 0.5, Qmin = 10.) 



(A = 0.5, Qmin = 10.) 




T(u)*H(u) 




0.4 0.5 0.6 0.7 



0.4 0.5 0.6 0.7 



(A = 0.5, Qmin = 20.0) 



(A = 0.5, Qmin = 20.) 




_ . . , n_> 

0.1 0.2 0.3 0.4 0.5 



(A = 0.5, independent of Qmin) 
Figure 2: The same as in the previous figure for A = 0.5 
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T(u) *H(u) 
0.7p 
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(A = 1.0, Qmin = 10.) 





0.4 0.5 0.6 0.7 



0.4 0.5 0.6 0.7 



(A = 1.0, Qmin = 20.0) 



(A = 1.0, Qmin = 20.) 




u-> 

0.1 0.2 0.3 0.4 0.5 



(A = 1.0, independent of Qmin) 
Figure 3: The same as in the previous figure for A = 1.0 



Table 3: The same quanitities with table I for A = 0.5 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t 


0.0 


1.690 


1.241 


0.861 


0.558 


0.453 


0.372 


0.224 


h 


0.0 


0.198 


0.151 


0.107 


0.083 


0.076 


0.071 


0.063 


t 


10. 


0.000 


0.267 


0.337 


0.268 


0.229 


0.194 


0.122 


h 


10. 


0.000 


0.344 


0.175 


0.113 


0.097 


0.087 


0.072 


h 


20. 


0.000 


0.000 


0.121 


0.123 


0.111 


0.098 


0.066 


h 


20. 


0.000 


0.000 


0.267 


0.150 


0.124 


0.106 


0.081 



Table 4: The same quanitities with table I for A = 1.0 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t 


0.0 


1.729 


1.299 


0.919 


0.600 


0.487 


0.399 


0.240 


h 


0.0 


0.314 


0.247 


0.181 


0.141 


0.131 


0.123 


0.112 


t 


10. 


0.000 


0.252 


0.353 


0.289 


0.247 


0.209 


0.132 


h 


10. 


0.000 


0.579 


0.291 


0.187 


0.163 


0.147 


0.124 


t 


20. 


0.000 


0.000 


0.120 


0.131 


0.120 


0.106 


0.071 


h 


20. 


0.000 


0.000 


0.455 


0.249 


0.205 


0.177 


0.137 



shown in Fig. 4 for A = 0, 1. We see that in the case of A = 1.0 the positive 
section of the function is enhanced. 

6.2 The Directional Rates 

Once again we distinguish two cases, the total and the differential rates. 
6.2.1 The Directional Total Event Rates 

The directional total event rates, which arize by summing the directional 
rates in all three directions is beyond the goals of the present experiments. 
We will, however, include it in the present discussion. The unmodulated rates 
can be can be parameterized in terms of the parameter t'^. This describes 
the modification of the total directional non modulated event rate due to the 
convolution with the velocity distribution. The modulation is now described 
by the three parameters /ii, /12, /i-3- These are shown in tables 5-7. It is clear 
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( /(y) ^Fi,Gi, Fi + Gi) ( /(y) ^F2,G2, A + G2 ) 

Figure 4: The quantities f{y) described in the text associated with Fi{X, 2(A + l)y), Gi(A, y) 
and -Fi{A, 2(A + l)y) + (^^(A, y) , i = 1,2. Thick solid line corresponds to i = 1, 2,the 
finest line to j = 1, 2 and the dashed line to the sum of the two. The intemediate thickness 
line corresponds to A = 0, in which case C?i = 0, i = 1, 2 



from Eq. ( |112| ) that the modulation of the total rate is no longer given by a 
simple sinusoidal function. For some interesting cases the situation is shown 
in Fig. 5. 

An idea about what is happening can be given by hm and am- The first 
gives the difference between the maximum and the minimum values of modu- 
lated amplitude. The second involves the phase shift from the second of June, 
which is no longer the the date of the maximum. 

The second of June gives the location of the maximum, when only the 
component of the earth's velocity along the sun's direction of motion is con- 
sidered or when /13 is neglected. In almost all cases considered in this work, 
however, /13 is important and in fact the obtained shift is on the average about 
± 35 days from the second of June. 

From tables 5-7 we see that without detector energy cutoff, Qmin = 0, 
t° is a decreasing value of the LSP mass. It takes values of about 2 for low 
LSP mass and is decreased by an order of magnitude as we go to higher LSP 
masses. This is due to the nuclear form factor effects, which are present but 
not so severe for this intermediate mass nucleus. It is not a sensitive function 
of the asymmetry parameter A. 

As expected, in the presence of energy cut off, is greatly reduced and 
becomes unobservable for light LSP. As the LSP mass increases increases, 
(see tables 5-7. It reaches a maximum at about 80 GeV and it starts decreasing. 
But even for the heaviest LSP this reduction is not much larger than 1/3, even 
for Qmin — 20KeV. In other words a heavy LSP can cause sufficient energy 
transfer to partially compensate for the loss of phase space. 

We also see that in all cases hm is much larger than 2hi, suggesting that 
when it comes to directional detection all components of the earth's velocity are 
important. We also notice that for those LSP masses, which give a detectable 
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Phase of Earth a — > rad 
(a) 



Phase of Earth a — > rad 
(b) 



Figure 5: The the total directional modulation amplitude as a function of the phase of 
the earth a in two cases. a) for hi = 0.059, /12 = 0.117 and hs = 0.135 and b) hi = 0.192, 
/i2 = 0.146 and /13 = 0.231. Note that in case b the minimum is negative. The results shown 
are for the target 53!^^''^ (for the definitions see text). 



total rate, the modulation amplitude does not appreciably change with the 
LSP mass. It is not greatly affected by the energy cutoff Qmm- It is enhanced, 
however, by about a factor of two in going from A = 0, no asymmetry, to A = 1, 
maximum asymmetry. 

From the above discussion it is clear that one needs all three modulation 
parameters, /ii,/i2,/i3. We remind the reader that in Eq (110) the z-axis has 
been chosen in the direction of the sun's velocity, the y-axis is perpendicular to 
the plane of the galaxy and the x-axis in the radial (galactocentric) direction. 
We mention again that /12 and /is are constant, 0.117 and 0.135 respectively, 
in the symmetric case. On the other hand hi and /13 substantially increase in 
the presence of asymmetry. 

The precise value of the directional rate depends on the direction of obser- 
vation. One can find optimal orientations, but we are not going to elaborate 
further. 



6.2.2 The Directional Differential Event Rates 



The directional differential rate is also very hard to detect, but perhaps a 
bit more practical than the total rates described in the previous subsection. It 
can be in terms of four functions of u, namely Rq{u) and Hi{u), i = 1,2, 3. 

The situation is rather complicated and following our discussion of the 
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Table 5: The quantities t",/ti and h,„ for A = in the case of the target 53/^^'' for vari- 
ous LSP masses and Qmin in KeV (for definitions see text). Only the scalar contribution 
is considered. Note that in this case /i2 and /13 are constants equal to 0.117 and 0.135 
respectively. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 




0.0 


1.960 


1.355 


0.886 


0.552 


0.442 


0.360 


0.212 


hi 


0.0 


0.059 


0.048 


0.037 


0.029 


0.027 


0.025 


0.023 




0.0 


0.164 


0.144 


0.124 


0.111 


0.107 


0.104 


0.100 


t° 


10. 


0.000 


0.365 


0.383 


0.280 


0.233 


0.194 


0.119 


hi 


10. 


0.000 


0.086 


0.054 


0.038 


0.033 


0.030 


0.025 




10. 


0.000 


0.214 


0.155 


0.127 


0.119 


0.113 


0.104 


t° 


20. 


0.000 


0.080 


0.153 


0.136 


0.11 


0.102 


0.065 


hi 


20. 


0.000 


0.123 


0.073 


0.048 


0.041 


0.036 


0.028 


hm 


20. 


0.000 


0.282 


0.190 


0.145 


0.132 


0.123 


0.109 



previous section we will give gross description of the modulation using the 
functions aij(u) and Hm{u). The phase shiftaij has been found to be a con- 
stant and about 0.7, which corresponds to a shift about ±35 days from June 
the second. Since , however, Hm is defined as the ratio of two quantities, it can 
appear large because the denominator (non modulated rate) becomes small. 
Thus, following the strategy of the previous subsection we also present the 
quantity Rm — RaH„i- These functions are shown in Fig. 6 for LSP masses in 
the range 30-250 GeV, A = and Qmin = 0, lOKeV and 20KeV. Note that the 
quantity Hm. is itself independent of the cutoff except that one should look at 
the u relevant to the allowed energy transfer interval. 

The curves shown correspond to LSP masses as in the undirectional case. 
Again due to the normalization of RO, the area under the corresponding curve 
is unity. 

The above quantities for A = 0.5 and 1.0 are shown in Fig. 7, but only 
for Qmin=0. Their dependence on the energy transfer cutoff shows behavior 
similar to that of Fig. 1. In any case for Qmin=^0, 20 KeV, the functions Rq 
and Rm show a behavior similar to that for Qmin = 0, except that they start 
from higher energy transfer. We should remind the reader, however, that that 
in all cases Rq represents the relative differential rate, i.e. it is normalized so 
that the area under the corresponding curve is unity for all LSP masses. One, 
therefore, should take into account the factor of tables 5-7. 
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Figure 6: The relative difTerential event rate RO and the amplitudes for modulation Rm and 
Hm vs u for the target 53/^^^ in the case of symmetric velocity distribution, A = (for the 
definitions see text). 
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Figure 7: The same as in Fig.2 in the asymmetric case (A = 0.5 and A = 1.0). Only the case 
Qmin = is exhibited. 
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Table 6: The same as in the previous table, but for the value of the asymmetry parameter 
A = 0.5. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


2.309 


1.682 


1.153 


0.737 


0.595 


0.485 


0.288 


hi 


0.0 


0.138 


0.128 


0.117 


0.108 


0.105 


0.103 


0.100 


h2 


0.0 


0.139 


0.137 


0.135 


0.133 


0.133 


0.133 


0.132 


ha 


0.0 


0.175 


0.171 


0.167 


0.165 


0.163 


0.162 


0.162 


hm 


0.0 


0.327 


0.307 


0.284 


0.266 


0.261 


0.257 


0.250 


tO 


10. 


0.000 


0.376 


0.468 


0.365 


0.308 


0.259 


0.160 


hi 


10. 


0.000 


0.174 


0.139 


0.120 


0.114 


0.110 


0.103 


h2 


10. 


0.000 


0.145 


0.138 


0.135 


0.134 


0.134 


0.133 


h3 


10. 


0.000 


0.188 


0.174 


0.167 


0.165 


0.164 


0.162 


hjji 


10. 


0.000 


0.400 


0.328 


0.290 


0.278 


0.270 


0.256 


tO 


20. 


0.000 


0.063 


0.170 


0.171 


0.153 


0.134 


0.087 


hi 


20. 


0.000 


0.216 


0.162 


0.133 


0.124 


0.118 


0.107 


h2 


20. 


0.000 


0.155 


0.143 


0.137 


0.136 


0.135 


0.133 


hs 


20. 


0.000 


0.209 


0.182 


0.171 


0.168 


0.166 


0.164 


hm 


20. 


0.000 


0.487 


0.374 


0.316 


0.299 


0.286 


0.265 



From the plots shown one can see that Rm can be quite large, about 20 % 
of Rq, but it falls slower as a function of the energy transfer. For this reason 
the modulation amplitude Hm is increasing as a function of u. It is interesting 
to note that the modulation amplitude Hm is increased by more than a factor 
of two, as the asymmetry parameter A changes from zero to one, for all energy 
transfers. Thus, even at zero energy transfer, for A = 1 the variation in the 
amplitude due to the earth's motion can increase by about 40% between the 
minimum (around December 2) and the maximum (around July 10 or the end 
of may), a big effect indeed. It can become even larger if one can restrict 
oneself to only part of the phase space, i.e. if one is satisfied with fewer counts. 

In the case of the directional differential rate one clearly needs, in addition 
to i?o , the functions Hi{u), 1=1,2,3, which are plotted in Fig. 4. In the case 
of A = only Hi is plotted, since the other two are in this case constant 
{H2 = 0.117, Hs = 0.135). We see that in the presence of asymmetry, e.g. 
A = 0.5 and 1.0, all functions, but especially Hi and H3, are substantially 
increased. 
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Figure 8: The same as in Fig.3 for the quantities H\,H2 and Hz. These quantities do not 
depend on Q^i„, except for the fact one should look at u > Umin- 
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Table 7: The same as in the previous, but for the value of the asymmetry parameter A = 1.0. 











LSP 




in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


tO 


0.0 


2.429 


1.825 


1.290 


0.837 


0.678 


0.554 


0.330 


hi 


0.0 


0.192 


0.182 


0.170 


0.159 


0.156 


0.154 


0.150 


h2 


0.0 


0.146 


0.144 


0.141 


0.139 


0.139 


0.138 


0.138 


h3 


0.0 


0.232 


0.222 


0.211 


0.204 


0.202 


0.200 


0.198 




0.0 


0.456 


0.432 


0.404 


0.382 


0.375 


0.379 


0.361 


tO 


10. 


0.000 


0.354 


0.502 


0.410 


0.349 


0.295 


0.184 


hi 


10. 


0.000 


0.241 


0.197 


0.174 


0.167 


0.162 


0.154 


hi 


10. 


0.000 


0.157 


0.146 


0.142 


0.140 


0.140 


0.138 


h3 


10. 


0.000 


0.273 


0.231 


0.213 


0.208 


0.205 


0.200 




10. 


0.000 


0.565 


0.464 


0.413 


0.398 


0.387 


0.370 


tO 


20. 


0.000 


0.047 


0.169 


0.186 


0.170 


0.150 


0.100 


hi 


20. 


0.000 


0.297 


0.226 


0.190 


0.179 


0.172 


0.159 


h2 


20. 


0.000 


0.177 


0.153 


0.144 


0.142 


0.141 


0.139 


h3 


20. 


0.000 


0.349 


0.256 


0.224 


0.216 


0.211 


0.203 




20. 


0.000 


0.709 


0.550 


0.448 


0.424 


0.408 


0.380 



7 Conclusions 

In the present paper we have expanded the the results obtained in of our 
recent letter c3. We have calculated all the parameters, which can describe 
the modulation of the direct detection rate for supersymmetric dark matter. 
The differential as well as the total event rates were obtained both for the 
non directional as well as directional experiments. All components of the 
earth's velocity were taken into account, not just its component along the 
sun's direction of motion. Realistic axially symmetric velocity distributions, 
with enhanced dispersion in the galactocentric direction, were considered. The 
obtained results were compared to the up to now employed Maxwell-Boltzmann 
distribution. 

We presented our results in a suitable fashion so that they do not depend on 
parameters of supersymmetry other than the the LSP mass. Strickly speaking 
the obtained results describe the coherent process in the case of ^^^/, but we 
do not expect large changes, if the axial current is considered. Recall that the 
dependence on supersymmetry is contained in the parameter R not discussed 
in the present paper. The nuclear form factor was taken into account and the 
effects of the detector energy cut off were also considered. 
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Our results, in particular the parameters t and to (see tables 2-4 and 5- 
8) indicate that for large reduced mass, the advantage of fir (see Eqs. (p^)- 
( [Iq ) is lost when the nuclear form factor and the convolution with the velocity 
distribution are taken into account. 

In the case of the undirectional total event rates we find that in the sym- 
metric case the modulation amplitude for zero energy cutoff is less than 0.07. 
It gets substantially increased in the case of asymmetric velocity distribution 
with largest asymmetry (A = 1). It can reach values up tO .31. In the presence 
of the detector energy cutoff it can increase even further up to 0.46, but this oc- 
curs at the expense of the total number of counts. The modulation amplitude 
in the case of the differential rate is shifted by the asymmetry at higher energy 
transfers and, for maximum asymmetry A = 1, gets about doubled compared 
to the symmetric case (A — 0). This amplitude does not depend on the the 
energy cutoff, but the lower energy transfers, will, of course, be excluded if 
such a cutoff exists. 

Analogous conclusions can be drawn about the directional differential 
event rate. The presence of asymmetry more than triples the differential mod- 
ulation amplitude (from about 10% to about 35%). There exist now regions of 
the energy transfer such that the modulation amplitude can become as large 
as 50%. 

Finally it is important that one should consider all components of the 
earth's motion, not just its velocity along the sun's motion, especially if the 
directional signals are to be measured. 
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